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ABSTRACT 


A  simple  expression  is  derived  for  the  influence  of  roll  acceleration  on 
the  angle -of-attack  convergence  of  rolling  reentry  vehicles  with  pitch  or  yaw 
damping.  Also  included  is  an  analysis  of  the  windward -meridian  rotation  rate 
of  the  rolling  vehicle,  which  is  coupled  through  the  gyroscopic  equations  of 
motion  to  the  angle  of  attack  and  roll  rate.  It  is  found  that  two  modes  of 
motion  exist  for  which  the  windward  meridian  is  generally  oscillatory  in  one 
mode  and  rotary  in  the  other  (i.  e.,  the  vehicle  either  oscillates  about  its  roll 
axis  relative  to  the  wind  or  rotates  continuously  in  one  direction).  The  mode 
of  motion  depends  on  the  initial  reentry  conditions,  and  it  is  shown  that  a 
roll-induced  instability  can  occur  whereby  the  motion  changes  from  the  oscil¬ 
latory  mode  to  the  more  stable  rotary  mode.  The  analytical  approximations 
are  compared  with  computer  solutions  of  the  complete  equations  of  motion. 
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1.  INTRODUCTION 


Of  particular  interest  in  the  performance  evaluation  of  a  reentry  vehicle 
is  its  angle-of-attack  behavior  during  reentry  through  the  atmosphere.  Under¬ 
standing  of  the  angle-of-attack  behavior  is  essential  for  prediction  of  the 
vehicle's  trajectory  and  for  assessment  of  the  thermal  requirements  of  the 
vehicle's  heat-protective  system.  In  addition  to  the  heating  prediction,  a 
knowledge  of  the  windward-meridian  rotation  rate  is  required  if  one  is  to 
couple  the  heating  or  ablation  rate  prediction  with  the  vehicle  dynamic  motion. 

This  paper  treats  the  influence  of  roll  rate  and  roll  acceleration  on  the 
angle-of-attack  convergence  of  rolling  reentry  vehicles  and  presents  a  tech¬ 
nique  for  predicting  the  windward -meridian  rotation  rate,  which  is  coupled 
through  the  gyroscopic  equations  of  motion  to  the  angle  of  attack  and  roll  rate. 
The  windward-meridian  rotation  rate  is  also  dependent  on  the  initial 
(exoatmospheric)  motion,  and  it  is  found  that  two  modes  of  motion  exist  for 
which  the  windward -meridian  is  generally  oscillatory  in  one  mode  and  rotary 
in  the  other  (i.  e. ,  the  vehicle  either  oscillates  about  its  roll  axis  relative  to 
the  wind  or  rotates  continuously  in  one  direction).  It  is  further  shown  that  the 
rotary  mode  is  more  stable,  although  either  mode  can  persist,  depending  on 
the  initial  conditions,  and  an  instability  can  occur  with  the  presence  of  roll 
acceleration  whereby  the  motion  changes  from  the  oscillatory  mode  to  the 
more  stable  rotary  mode.  A  criterion  is  presented  for  predicting  this  insta¬ 
bility,  which  occurs  only  when  the  roll  rate  is  large  relative  to  the  natural 
pitch  frequency  of  the  vehicle. 

The  angle-of-attack  convergence  behavior  of  rolling  reentry  vehicles 
has  been  studied,  and  a  number  of  papers  on  the  subject  are  reported  in  the 
literature.  *  All  of  these,  however,  treat  the  case  of  constant  roll  rate; 
only  one  includes  the  effect  of  pitch  or  yaw  damping  on  the  convergence 
envelope,  *  and  none  discusses  in  any  detail  the  influence  of  exoatmospheric 
conditions  on  the  subsequent  angle-of-attack  convergence  and  vehicle  motion. 
Because  of  the  susceptibility  of  high-performance  reentry  vehicles  to  roll 
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resonance  and  other  phenomena  that  can  cause  significant  roll  rate  excursions, 
it  Is  of  interest  to  assess  the  influence  of  roll  acceleration  on  the  angle -of- 
attack  convergence  and  on  the  windward -meridian  rotation  behavior. 

The  present  study  treats  the  vehicle  motion  during  the  period  in  which, 
because  of  misalignment  with  the  flight  path  during  reentry,  the  angle  of 
attack  is  large  relative  to  the  trim  angle  of  attack  from  mass  or  configura¬ 
tional  asymmetries.  The  quasi-steady  or  mean  value  of  the  angle -of-attack 
oscillation  envelope  is  obtained  as  a  function  of  time  (or  altitude)  during  this 
period. 


II.  EQUATIONS  OF  MOTION 


The  vehicle  motion  is  described  in  terms  of  the  Euler  angles  >&,  0, 

as  shown  in  Fig.  1.  If  the  principal  moments  of  inertia  about  the  %,  q,  { 
axes  are  I,  I,  Ix,  respectively  (pitch,  yaw,  roll),  and  the  aerodynamic 
moments  about  these  axes  are  Mt,  M  ,  M_,  the  moment  equations  of  motion 
in  terms  of  the  Euler  angles  may  be  written 


M*  =  10  +  i  pi  sin  6  -  ii  sin©  cos  0 
9  xr 

M  =  I  4:  (i  sin  0)  +  10*  cos  0  -  I  p0 
T]  dt  x 

M,  =  1 
L  x  dt 


(1) 


where  the  roll  rate  p  is  defined  by 


p  =  i  +  i  cos  0, 


(2) 


The  angular  rate  4>  is  tne  windward-meridian  rotation  rate,  i.  e. ,  the  roll  rate 
with  respect  to  the  wind.  It  is  assumed  that  the  aerodynamic  moments  consist 
only  of  pitch  or  yaw  moments  from  angle  of  attack,  an  arbitrary  roll  moment, 
and  pitch  and  yaw  damping  moments.  It  is  further  assumed  that  the  vehicle  is 


acrodynamically  axisymmetric  (i.  e.,  CM  =C 

IN/v 


„  ),  so  that  the  restoring  torque 
from  angle  of  attack  is  independent  of  the  roll  orientation  4>  and  dependent  only 
on  0.  The  moments  are  then  written 
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*  sin  8 


=  C^Sd 


(3) 


where  f  is  the  vehicle  static  margin  (considered  positive  for  the  statically 
stable  vehicle),  and  the  damping  terms  consist  of  both  the  pitch  (or  yaw) 
damping  derivative  Cm^  plus  a  normal  force  damping  term  2C^  I/md 
to  account  for  lateral  motion  of  the  vehicle  center  of  mass.41 

Substituting  Eq.  (3)  in  Eq.  (i)  and  dividing  the  first  two  equations  by  I 
and  the  third  by  Ix  gives 

2  .  •  •  2 
-u>  8  -  v0  =  6  +  tipf  sin  9  -  *  sin  6  cos  6 

-vtf  sin  0  =  ^  (*  siu  6)  +  0*  cos  6  -  np0  (4) 

*  =  p 

where  the  coefficients  are 


The  parameter  u>  is  the  natural  pitch  frequency  of  the  vehicle. 


$ 

It  can  be  shown,  for  an  axisymmetric  vehicle,  that  the  lateral  motions  of  the 
center  of  mass  from  normal  forces  acting  at  the  center  of  pressure  are  pro* 
portional  to  pitch  or  yaw  rotations  about  a  point  a  distance  I /ml  ahead  of  the 
center  of  mass. 
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Ul.  QUASI -STEADY  SOLUTION  FOR  ANGLE -OF -ATTACK 

CONVERGENCE 


The  vehicle  motion  during  reentry,  described  in  terms  of  the  time- 
dependent  Euler  angles  [  4(t),  ^(t),  0 (t)]  ,  is  assumed  to  be  of  the  form 

[*<t).*(t)fe(t)]  =  [*(t),f(t),0(t)]  +  [*+(t),*+(t),e+(t)]  (6) 

where  [  *(t),  f(t),7(t)]  represents  a  quasi-steady  component  that  varies  rela¬ 
tively  slowly  with  time  (of  the  order  of  the  dynamic  pressure)  and 
[¥+(t),  ^+(t),  0+(t)]  represents  an  oscillation  of  higher  frequency  about  the 
average  (quasi-steady)  values.  In  the  ensuing  analysis,  the  oscillations  are 
neglected  and  a  solution  is  obtained  for  the  quasi-steady  angle -of-attack  con¬ 
vergence  behavior  of  the  descending  vehicle.  This  development  is  heuristic 
in  nature  and  based  on  physical  intuition.  It  is  justified  on  the  basis  of  agree¬ 
ment  with  more  complete  and  lengthier  analyses  presented  later. 

Taking  sin  0*0  and  cos  0*1  and  neglecting  pitch  damping  v  in  the  first 
of  Eqs.  (4)  gives,  for  the  pitch  equation, 

0  +  (w2  +  pptf  -  *2)  e  =  o  (7) 

This  equation  describes  a  nonlinear  oscillation  in  0  about  a  nonzero  quasi¬ 
steady  value  9,  since  0,  by  definition,  is  always  positive.  One  can  interpret 

T  „ 

the  corresponding  quasi-steady  values  of  ¥to  be  those  that  make  0  in  Eq.  (7) 
zero,  thus  satisfying  the  relation 

2  **■  t-2 

u  +  pp¥  -  =0 

or 


*  =  (M-p/2)  ±  [  (n-p/2)2  +  u2]  *^2 


(8) 


This  result  is  valid  as  long  as  9  is  appreciably  greater  than  the  trim  angle  of 
attack  from  aerodynamic  or  mass  asymmetries  (assumed  to  be  zero  for  the 
present  analysis).  For  y.p/2  small  relative  to  u>,  the  two  values  for  may  be 
'V  ‘tten  approximately 


(fr  •-  w  +  -  w  +  (9) 

This  indicates  that,  for  values  of  angle  of  attack  sufficiently  greater  than 
trim,  the  precession  rate  can  oscillate  in  either  a  positive  or  negative  mode, 

the  quasi-steady  value  of  which  is  slightly  displaced  from  the  natural  pitch 

£ 

frequency  for  a  slender  vehicle,  as  shown  in  Fig.  2.  The  mode  of  oscillation 
depends  on  the  initial  reentry  conditions  and  is  discussed  later.  Once  i  is 
determined,  the  windward-meridian  rotation  for  small  0,  is  simply  the 
difference  between  the  roll  rate  p  (assumed  known)  and 

The  angle -of-attack  convergence  can  be  obtained  from  the  second  of 
Eqs.  (4).  This  equation  can  be  written 

('i'  sin20)  +  vi  sin20  -  pp9  sin  6=0  (10) 

which,  for  small  9,  can  be  differentiated  to  give 

(i-‘?)T,(E’2)+(#+‘"i')92  =  0-  <"> 


Th*  two  precession  modes,  Eqs.  (8)  or  (9),  are  analogous  to  the  fast  and 
slow  precession  of  a  spinning  top  but  differ  in  sign,  because  the  pitch  moment 
acts  in  the  opposite  direction  to  the  top  turning  moment  for  a  statically  stable 
missile. 


Wc  now  consider  the  time  variation  of  the  quasi-steady  values  0  and  ¥  by 
writing  Eq.  (11)  in  the  form 


^-e2  +  f(t)62  =  0, 

where 

f(t)  s 

Equation  (12)  has  the  familiar  solution 


(dj'/dt)  +  v'j' 
*  -  C^p/  2) 


(12) 


(13) 


(14) 


On  substituting  ^  from  Eq.  (8)  in  Eq.  (13)  and  integrating  this  result  in 
Eq.  (14)  we  find,  for  the  quasi-steady  angle-of-attack  convergence,  the 
expression 

1  -t  .  f  ?  -1/2]  | 

-J  J  (£  +  vjj^l  ±  (1  +  O  Jdt  (15) 

where  o- =  2u/iip.  The  plus  and  minus  signs  in  the  exponent  correspond  to  the 

positive  and  negative  modes,  respectively,  of  the  precession  rate  4r,  depicted 

in  Fig.  2.  The  roll  acceleration  term  p/p  in  the  exponent  is  an  effective 

damping  term  that  adds  to  or  subtracts  from  the  damping  coefficient  v,  depend 

£  2  .1/2 

ing  on  the  sense  of  the  roll  acceleration.  Similarly,  the  term  (1  +  <r  } 


t  =  <•  +*2> 

o 


It  appears  that  negative  (decelerating)  values  of  p  that  yield  negative  values 
for  the  effective  damping  coefficient  [  (p/p)  +  v]  would  produce  a  momentary 
angle-of-attack  divergence.  Such  cases  might  prove  interesting  for  future 
investigations . 


either  increases  or  decreases  the  effective  damping  and,  therefore,  increases 

or  decreases  the  rate  of  angle -of -attack  convergence,  depending  on  whether 

the  precession  mode  is  positive  or  negative.  The  conditions  that  determine 

the  precession  mode  are  discussed  later. 

For  a  slender  reentry  vehicle.  |A <c  i ,  the  parameter  cr^  is  generally 

much  greater  than  unity  except  for  highly  supercritical  roll  rates  p  »u. 

2  -1/2 

Therefore,  in  most  cases  the  term  ( i  +  or  )  in  the  exponent  of  Eq.  (15) 

can  be  ignored  compared  with  unity,  and  Eq.  (15), with  the  assumption 
2  2 

1  +  <r  *  or  ,  reduces  to 


(16) 


For  an  exponential  atmosphere  of  scale  height  H,  and  with  a  straight-line 
trajectory  of  path  angle  y,  the  integral  of  Eq.  (16)  can  be  evaluated  and  the 
angle -of -attack  convergence  ratio  reduces  to  the  simple  result 


where 


(i?) 


(i8) 


and  p  is  density.  The  results,  Eq.  (16)  or  (17),  indicate  that  for  roll  rates 
that  are  not  excessively  supercritical  [i.  e. ,  for  (i  +  <r  )"  '  «i)] ,  the  angle- 

of -attack  convergence  depends  only  on  the  reentry  roll  rate  and  is  independent 
of  roll  acceleration. 

For  the  constant  roll  rate  case  (p  =  o),  the  result,  Eq.  (15),  reduces  to 
an  expression  similar  to  that  derived  in  Ref.  i  for  the  negative  precession 
mode. 


IV.  INITIAL  REENTRY  CONDITIONS 


Before  reentry,  the  vehicle  is  in  a  state  of  moment -free  motion,  and 
the  relation  between  the  various  angular  rates  is  uniquely  determined.  In 
general,  the  vehicle  will  have  some  initial  roll  rate  pQ  and  will  undergo  a 
steady  precession  (coning  motion)  with  a  precession  rate  ft  and  cone  half¬ 
angle  a  (Fig.  3).  The  cone  axis  will,  in  general,  be  inclined  to  the  flight  path, 
and  the  velocity  vector  may  lie  inside  or  outside  the  coning  circle.  The  latter 
case  is  shown  in  Fig.  3.  The  exoatmospheric  roll  rate,  precession  rate,  and 
cone  half-angle  are  related  by  the  expression 


cos  a 


(19) 


This  relation  follows  from  the  first  of  Eqs.  (4)  with  6  =  0  =  w  =  o,  and  a  and  ft 
substituted  for  e  and  tf,  respectively. 

After  entering  the  atmosphere,  the  statically  stable  vehicle  is  subjected 

to  an  aerodynamic  pitch  moment  that  tends  to  align  the  vehicle  with  the  flight 

path.  However,  because  of  the  angular  momentum  comprised  of  the  roll  and 

precession  rates,  which,  on  the  average,  is  directed  along  the  coning  axis, 

the  pitch  motion  is  resisted  by  gyroscopic  forces  that  induce  a  precession  of 

the  angular  momentum  vector  about  the  average  flight  path.  This  precession 

is  opposite  in  direction  to  the  angular  momentum  (retrograde  precession).  A 

projection  of  the  path  described  would  e  as  shown  in  Fig.  4  when  viewed  along 

the  direction  of  flight.  The  residual  coning  motion,  in  the  same  direction  as 

the  angular  momentum  vector  (e.g.,  clockwise),  has  been  called  nutation  by 
7 

Nicolaides,  and  the  retrograde  precession  in  the  opposite  direction  (counter - 
clockwise)  has  been  called,  simply,  precession.  In  general,  the  two  motions 
exist  simultaneously. 


These  definitions  are  more  restrictive  than  are  the  classical  definitions  of 
nutation  and  precession  as  being  variations  in  6  and  ¥,  respectively  (see,  for 
example,  Ref.  8,  p.  432). 
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The  expression  derived  above  for  angle -of-attack  convergence,  Eq.  (15), 
corresponds  to  two  limiting  cases  of  the  type  of  motion  depicted  in  Figs.  3 
aad  4.  The  first  case,  with  the  negative  sign  in  the  exponent,  represents  the 
limiting  condition  of  aero  coning  half -angle  for  which  the  vehicle  is  initially 
at  angle  of  attack  with  respect  to  the  flight  path  and  is  rolling,  as  shown  in 
Fig.  5a.  The  second  case,  with  the  positive  sign  in  the  exponent,  is  the 
condition  in  which  the  axis  of  the  initial  coning  motion  coincides  with  the 
average  flight  path,  as  shown  in  Fig.  5b. 

The  Euler  angles  (ty  f ,  0)  were  defined  with  respect  to  am  inertial  frame 
of  reference  that  moves  along  the  flight  path  such  that  i  represents  precession 
about  the  velocity  vector  (Fig.  1).  Therefore,  the  two  exoatmospheric  condi¬ 
tions  represented  in  Figs.  5a  and  5b  correspond  to  initial  conditions  on  *  of 
4^  *  0  and  =  ft  *  ppo/cos  0O,  respectively.  Referring  back  to  the  quasi¬ 
steady  solution  for  ♦,  Eq.  (8),  we  note  that,  in  the  limit  as  u  —  0,  correspond¬ 
ing  to  aero  atmospheric  density,  $  can  have  the  two  values  0  and  p-p  (the  latter 
value  would  be  ^p/cos  6,  without  the  small  angle  approximations  in  Eq.  7). 
Therefore,  the  two  precession  modes  shown  in  Fig.  2  follow  from  the  two 
cases  of  exoatmospheric  motion  shown  in  Fig.  5.  Under  certain  conditions 
when  the  initial  precession  mode  is  positive,  an  instability  can  occur  in  the 
presence  of  roll  acceleration  whereby  the  motion  changes  from  the  positive  to 
the  negative  precession  mode.  A  criterion  for  predicting  this  instability  is 
derived  in  the  next  section. 
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V.  PRECESSION  INSTABILITY 


Equation  (11)  describes  the  relation  between  the  total  angle  of  attack  0 
and  the  precession  rate  4%  without  any  assumption  of  quasi-steady  motion. 

This  equation  can  be  written  in  the  form 

0  =  .  (d*/dt)  +  (20) 

2$  -  pp 

which  reveals  much  about  the  vehicle  motion.  First,  we  note  that  8,  by  defini¬ 
tion,  is  always  positive  and  4',  from  the  above  considerations,  oscillates  about 
either  the  positive  or  the  negative  branch  of  the  natural  pitch  frequency 
curve  co,  as  shown  in  Fig.  2.  Therefore,  for  4  >pp/2,  the  denominator  of 
Eq.  (20)  remains  positive,  and,  for  small  damping,  d6/dt  =  0  at  approximately 
d4/dfc  =  0,  so  that  6  and  4  oscillate  at  the  same  frequency  with  the  6  maxima 
corresponding  to  the  4  minima  and  vice  versa.  Similarly,  for  4  <  H-p/2,  the 
maxima  and  minima  of  both  6  and  4  are  in  phase,  as  shown  in  Fig.  6.  These 
results  can  be  deduced  in  another  manner.  Equation  (10),  with  the  assump¬ 
tions  v  =  0  and  sin  0  ~  8,  can  be  written  in  the  form 

a1,'*®2’-1?;^62^0  <21> 

With  the  further  assumption  that  p  is  constant  (or  slowly  varying  relative  to 
the  6  or  4  oscillations),  Eq.  (21)  can  be  integrated  to  give 

(4  -  ^)  0^  =  constant  (22) 


Equation  (22),  being  the  first  integral  of  an  equation  of  motion,  is  an  energy 
expression  and  indicates  that  |4  -  (H-p/2)  |  is  maximum  for  6  minimum  and 
vice  versa,  as  depicted  in  Fig.  6.  However,  if  4  is  in  the  positive  mode  such 
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Figure  6.  Angle -of - 

Oscillatioi 


that  ¥  >  H-p/ 2,  and  p  increases  because  o f  an  applied  roll  torque,  an  instability 
will  occur  as  h-p /2  approaches  the  minimum  value  of  the  4r  oscillation,  which 
causes  the  denominator  of  Eq.  (20)  to  approach  zero.  It  is  shown  from  a 
phase  plane  analysis,  which  follows,  that  at  the  point  of  instability  the  * 
oscillation  diverges  without  limit  and  finally  reverses  sign  from  an  infinitely 
large  positive  upper  bound  to  an  infinitely  large  negative  lower  bound.  There¬ 
after,  i  remains  in  the  more  stable  negative  mode.  The  divergence  in  4r 
corresponds  to  a  lower  bound  on  the  8  oscillation  that  approaches  zero,  from 
Eq.  (22). 

Equations  (7)  and  (20)  can  be  combined  to  give  a  single  equation  in  'fr 


(*  •  +  v  lr +  *  3(ir)  +  (,Ap  *  4v*>  ir + 


-  ^  J  +  h-p*  - 


(23) 


We  can  put  this  equation  in  a  form  suitable  for  phase  plane  analysis  by  making 
the  substitutions 


y  = 


.  dy  n  • 

V=y'dT"2P 


(24) 
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in  which  $  has  bstn  assumed  to  be  negligible.  Equation  (23)  is  then  written 


-frr2<y2  - 12)  +  cy(y  +  <*)  +(3/2Mv  +  b)  +  wy, 

dy  yv 


(25) 


where 

a*  [«2  +  (|ip/2)2]1/2 
b  s(2/3)|i  (*  +  vp) 

c  «d/2)(v2  -  2v) 


d«KP 


(26) 


The  parameters  a,  b,  c,  d,  and  v  vary  slowly  with  time  (of  the  order  of  the 
dynamic  pressure)  relative  to  y  and  are  assumed  to  be  constants  in  the  ensuing 
analysis.  The  stability  conditions  so  obtained  from  a  phase  plane  evaluation 
of  Eq.  (25)  apply  at  any  instant  of  time  with  the  appropriate  values  of  these 
parameters.  However,  since  these  parameters  do  change  with  time,  the 
nature  of  the  motion  and  the  stability  criteria  change  accordingly.  This  is 
discussed  in  greater  detail  later. 

Singularities  in  the  v  -  y  plane  (points  where  the  right  side  of  Eq.  (25) 
takes  the  form  sero  over  sero)  are  seen  to  occur  at 


y  =  0 
y  =  0 

y  =  yj 

y  =  y2 
y  =  y3 


V  =  0 

v  s  -b 
v  =  0 
v  =  0 
v  =  0 
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where  y^,  y 2>  and  y^  are  roots  of 

y3  -  (a2  +  f)y  -  T-  =  0.  (27) 

2 

The  parameters  c  and  d  are,  in  general,  much  smaller  than  a  ,  and  the  roots 
of  Eq.  (27)  can  be  shown  to  be  approximately  given  by 


a  a  +  *  + 


f  d 


y  »  -  a  -  «  + 


all +(«/»)] 
cd 


a[  1  +(« /a)] 


(28) 


-2f  d 

y3  ~ - y 

a[  l+(*/a)]  2 


where  <  is  small  relative  to  a  and  is  defined  by 


c 


(*2  “  2w) 

8[  w2  +  (np/2)2] 


m 


(29) 


We  now  examine  the  character  of  each  singularity  in  order  to  construct  the 
solution  curves  to  the  equation  of  motion.  Equation  (25)  can  be  written  in 
the  approximate  form  4 


dv 

dy 


-  2y  (y  -  a)  (y  +  a)  [  y  +  (2cd/a)]  +  j  v  (v  +  b)  +  vvy 

— 


(30) 


in  which  the  terms  containing  «  have  been  neglected  relative  to  a,  since  they 
do  not  change  the  character  of  the  singularity.  Consider  first  the  singularity 
at  y  =  a,  v  =  0.  If  the  origin  is  shifted  to  the  singularity  with  the  substitution 
x  =  y  -  a,  the  slope  dv/dx  in  the  limit  as  x,  v  -*0  is  found  to  be  of  the  form 
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Equation  (31)  is  a  special  case  of  the  differential  equation 


dv  _  Ax  +  Bv 


which  has  the  characteristic  equation 


(32) 


X2  -  X  (B  +  C)  -  (AD  -  BC)  =  0.  (33) 

The  nature  of  the  singularities  is  determined  by  the  relative  magnitudes  of 
the  constants  A*  B,  C,  and  D,  and  the  various  singularities  that  can  arise 
from  Eq*  (32)  have  been  previously  investigated.  Table  1  characterises 
these  singularities.  The  singularity  at  y  =  a,  v  =  0  described  by  Eq.  (31)  is 
characterised  by  the  coefficients 


Table  1.  Characteristics  of  the  Singularities  of  ^  + 


(B  -  C)2  +  4AD  >  0 
(B  -  C)2  +  4AD  <  0 

(B  -  C)2  +  4AD  =  0 


Node  if  AD  -  BC  <  0 
Saddle  if  AD  -  BC  >  0 


Center  if  B  +  C  *  0 
Spiral  if  B  +  C  *  0 


Node 


Stable  if  B  +  C  <  0 
Unstable  if  B  +  C  >  0 

Stable  if  B  +  C  <  0 
Unstable  If  B  +  C  >  0 

Stable  if  B  +  C  <  0 
I  Unstable  if  B  +  C  >  0 


A  =  -  4* 2  C  =  0 

B  =  ||  +  v  D  =  1 . 

in  which  the  constants  a,  b,  and  v  are  positive  by  definition.  For  this  case, 
(B  -  C)2  +  4AD  =  [  (3b/2a)  +  v]  2  -  16a2  <  0  (since  —  ,  v  «.  a)  and  B  +  C  = 

A 

(3b/ 2a)  +  v  >  0.  Therefore,  the  singularity  is  an  unstable  spiral.  Similarly, 
for  the  singularity  at  y  =  -a,  v  =  0,  if  the  origin  is  shifted  to  the  singularity 
with  the  substitution  x  =  y  +  a,  the  limiting  value  of  the  slope  dv/dx  is  found 
to  be 


lim  _ 

3[  0dv  -  -  4a  x  -  [(3b/2a)  -  v]v  (34) 

v-O^ 

For  this  case,  (B  -  C)2  +  4AD  =  [(3b/2a)  +  v]2  -  16a2  <  0,  as  before,  but 
B  +  C  =  -  (3b/2a)  +  v,  which  may  be  greater  or  less  than  zero,  depending  on 
the  damping  and  roll  acceleration.  Therefore,  this  singularity  is  either  a 
stable  or  unstable  spiral. 

The  quantity  B  +  C,  which  determines  the  degree  of  stability  or  insta¬ 
bility  of  the  spiral  singularities  at  ±a,  0  is  proportional  to  the  pitch  damping 
and  roll  acceleration,  which  are,  in  general,  small  quantities.  For  the  case 
of  zero  damping  and  constant  roll  rate  (i.  e. ,  v  =  p  =  0),  then  B  +  C  =  0  and 
the  spiral  singularities  reduce  to  center-type  singularities  in  this  limiting 
case.  The  singularities  at  ±a,  0  therefore  signify  the  existence  of  either 
stable  or  quasi-stable  oscillations  in  y  about  the  equilibrium  values  y  =  ±a, 
depending  on  the  magnitude  of  pitch  damping  and  roll  acceleration.  Moreover, 
the  oscillations  about  y  =  +a  are  of  either  constant  or  diverging  amplitude, 
corresponding  to  the  center  or  unstable  spiral  singularities,  whereas  the 
oscillations  about  y  =  -a  may  be  either  of  constant,  diverging,  or  converging 
amplitude,  depending  on  the  relative  magnitudes  of  the  pitch  damping  and  roll 
acceleration.  From  the  definition  of  y,  Eq.  (24),  the  oscillations  in  y  about 
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Aft  correspond  to  the  direct  and  retrograde  precession  oscillations  in  ¥  about 
the  positive  and  negative  branches,  respectively,  of  the  natural  pitch  frequency 
corves  identified  earlier  in  Fig.  2. 

The  above  considerations  are  based  on  the  assumption  of  quasi- steady 
motion  for  which  the  parameters  defined  by  Eq.  (26)  were  assumed  to  be 
constants.  However,  since  these  parameters  change  slowly  with  time,  they 
inflaeace  the  character  of  the  motion.  Consider,  for  example,  the  oscillation 
in  y  about  the  singularity  y  *  +a  and  let  there  be  a  positive  roll  acceleration 
so  that  the  roll  rate  p  increases  with  time.  The  increase  in  p  reduces  y, 
from  the  definition,  Eq.  (24),  whereas  the  value  a  about  which  y  oscillates 
increases  approximately  as  the  natural  pitch  frequency.  The  net  effect  is  to 
increase  the  amplitude  of  the  oscillation:  this  increase  is  characteristic  of 
the  unstable  spiral  singularity  that  exists  at  y  =  +a  with  the  presence  of  a  roll 
acceleration.  It  will  be  shown  that  if  p  increases  sufficiently  to  drive  y  to 
sero  (i.e. ,  pp/2  approaches  the  minimum  value  of  4).  an  instability  occurs 
whereby  the  motion  reverses  from  an  oscillation  about  y  =  +a  to  a  more  stable 
oscillation  about  y  =  -a.  The  nature  of  the  crossover  becomes  apparent  from 
the  character  of  the  remaining  singularities  at  and  near  the  origin. 

Consider  the  singularity  at  y  =  0,  v  =  -b.  The  origin  is  shifted  to  the 
singularity  by  the  substitution  w  =  v  +b.  The  slope  dw/dy  in  the  limit  as  y, 
w-*0  is 


lim 

dw 


t  v  .(4<ad/b)]  y  +  ^w 


/ 

(35) 


This  singularity  is  characterised  by  (B  -  C)2  +  4AD  -  1/4  >  0,  AD  -  BC 
*  -  3/2  <  0  and  B  +  C  *  5/Z  >  0,  which  is  an  unstable  node.  For  the  singu¬ 
larity  y  *  -  (2cd/a),  the  origin  is  shifted  to  the  singularity  with  the  substitu¬ 
tion  x*y  +  (2«d/a)  and  the  limiting  slope  is 


lim 


(36) 


x-0 

v-*0 


dv  _  2a^x  -  [(3ab/4id)  -  w]  v 
clx  "  v 


2  2  2 

The  characteristic  values  are  (B  -  C)  +  4AD  =  [  (3ab/4«d)  -  v)  +  8a  >0  and 
2 

AD  -  BC  =  2a  >0.  Therefore,  the  singularity  is  a  saddle.  Finally,  the 
singularity  at  the  origin,  y  =  0,  v  =  0,  has  the  limiting  slope 


lim 

dv_  4ad<y  +  (3/2)  bv  +  w  yv 
v^o  ** '  ** 


(37) 


This  is  a  higher-order  singularity  than  those  characterized  by  Eq.  (32),  and 
the  nature  of  the  singularity  is  not  readily  determined  analytically.  Conse- 
quently,  a  numerical  technique  such  as  the  method  of  isoclines*®  is  required 
for  determining  the  behavior  of  the  solution  curves  in  the  neighborhood  of  the 
singularity. 

Figure  7  summarizes,  schematically,  the  character  of  the  four  singu¬ 
larities  analyzed  above  plus  the  higher-order  singularity  at  the  origin  deduced 
from  limiting  cases  of  Eq.  (37).  The  relative  spacing  between  the  singularities 
is  grossly  distorted,  since  <  is,  in  general,  a  very  small  quantity  compared 
with  a.  The  qualitative  behavior  of  the  solution  curves  in  the  neighborhood  of 
the  singularities  near  the  origin,  as  determined  from  a  plot  of  Eq.  (25)  by  the 
method  of  isoclines,  is  shown  in  Fig.  8.  The  nature  of  the  crosoover  from  an 
oscillation  about  the  positive  branch  of  the  pitch  frequency  curve  to  an  oscil¬ 
lation  about  the  negative  branch  now  becomes  apparent.  Figure  8  shows  that 
there  are  two  switching  lines  in  the  positive  halfplane  and  one  in  the  negative 
halfplane.  To  the  right  of  the  first  switching  line,  an  oscillation  about  the 
positive  branch  persists  until  the  amplitude  reaches  sufficient  magnitude  to 
reach  the  first  switching  line  as  y  -*  0.  It  then  appears  possible  to  follow  a 
trajectory  between  the  two  switching  lines  which  will  pass  through  the  singu¬ 
larities  but  continue  around  the  positive  branch  (trajectories  labeled  3  and  4). 
Beyond  the  second  switching  line  the  trajectory  loops  around  the  negative 
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Figure  8.  Solution  Curves  in  Neighborhuod  of  Origin 
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branch  (y  =  -a)  and  can  return  to  the  positive  branch  provided  it  passes  to  the 
right  of  the  switching  line  in  the  negative  halfplane  (trajectories  1  and  2). 

This  motion  through  the  singularities  at  the  origin  represents  a  brief  period 
instability,  which  occurs  only  at  the  instant  y  —  0.  As  was  pointed  out 
earlier,  the  spacing  between  the  singularities  in  the  vicinity  of  the  origin  is 
extremely  small  compared  with  the  amplitude  of  the  motion.  Also,  in  view  of 
the  decreasing  value  of  y  due  to  the  increase  in  roll  rate,  it  would  be  expected 
that  only  a  small  number  of  trajectories,  if  any,  would  pass  through  this 
region.  To  the  left  of  the  origin,  as  y  continues  to  decrease  due  to  the 
increasing  roll  rate,  the  oscillations  appear  to  become  more  stable  about  the 
negative  branch,  as  represented  by  the  trajectories  to  the  left  of  the  singular¬ 
ity  at  -c .  This  behavior  has  been  verified  by  computer  solutions  of  the  equa¬ 
tions  of  motion,  which  are  discussed  later. 

Some  further  insight  can  be  obtained  from  a  physical  interpretation  of 
the  relation,  Eq.  (20),  between  the  angle  of  attack  and  the  precession  rate. 
From  the  definition  for  y,  Eq.  (24),  the  limiting  case  y  =  0  corresponds  to  a 
zero  value  of  the  denominator  on  the  right  side  of  Eq.  (20).  But  since  6,  by 
definition,  is  always  positive,  this  would  represent  an  instability  in  8  (i.  e. , 
a  large  value  of  6)  unless  either  8  or  (dv/dt)  +  vqr  also  approached  zero  along 
with  y.  What  is  found  from  the  computer  solutions  is  that  8  does  approach 
zero  in  the  region  of  instability,  since  the  6  oscillations  are  coupled  with 
the  '('oscillations  through  Eq.  (22),  as  depicted  in  Fig.  6.  From  conservation 
of  energy,  as  8  approaches  its  minimum  |4'  -  (p-p/ 2)  {  approaches  its  maximum 
and  vice  versa.  Therefore,  the  extreme  values  of  4'  (or  y),  both  positive  and 
negative,  will  occur  in  the  region  of  instability  as  y  -*0.  This  is  shown 
qualitatively  in  Fig.  9. 

The  foregoing  discussion  has  dealt  with  the  precession  rate  4f,  which  is 
related  to  the  roll  rate  p  and  the  windward-meridian  rotation  rate  p  through 
the  relation,  Eq.  (2), 


p  =  p  +  4  cos  8 
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(2) 


or.  for  siraii  A, 


p  *  p  +  4  .  (38) 

From  the  third  of  Eqs,  (4),  the  roll  rate  ia  gyros copically  uncoupled  from  the 
dther  angular  rates  and  depends  only  on  the  roil  moment  coefficient.  There¬ 
fore,  for  a  prescribed  roii  history,  the  windward-meridian  rotation  rate  is 
•imply  the  difference  p  =  p  -  4.  Depending  on  the  roll  history  and  on  whether 
4  is  in  the  positive  or  negative  mode,  the  windward  meridian  may  be  either 
"osciilatory"  or  "rotary";  i.e.  ,  the  vehicie  either  oscillates  about  its  roll 
axis  relative  to  the  wind  or  rotates  continuously  in  one  direction.  This  is 
illustrated  in  Fig.  10  for  an  arbitrary  roll  rate  history.  When  4  is  in  the 
positive  mode  with  the  roii  history  shown,  p  oscillates  from  positive  to 
negative  with  the  quasi-steady  value  $  •  p  -  u  -  «.  However,  for  4  in  the 
negative  mode,  4  does  not  change  sign,  and  the  windward  meridian  therefore 
rotates  continuously  in  one  direction.  Since,  as  was  shown  above,  the 
negative  4  mode  is  more  stable,  the  windward  meridian  is  more  likely  to  be 
rotary,  but  either  mode  may  persist,  depending  on  the  roll  history. 


VI.  COMPUTER  SOLUTIONS 


Figures  11  end  12  show  computer  solutions  of  the  equations  of  motiou, 
Eqs.  (4),  for  the  angle  of  attack  and  angular  rates  corresponding  to  the  two 
limiting  cases  of  exoatmospheric  motion  shown  in  Fig.  5.  Also  shown  for 
comparison  with  the  computer  results  are  quasi-steady  approximations  for 
the  angle-of-attack  convergence,  calculated  from  Eq.  (15).  Figure  11  corres¬ 
ponds  to  Fig.  5a,  in  which  the  rolling  vehicle  is  at  angle  of  attack  with  respect 
to  the  flight  path,  but  has  no  precession.  Figure  12  corresponds  to  the  other 
limiting  case,  in  which  the  vehicle  is  precessing  initially  about  the  average 
flight  path  in  the  positive  direction. 

Since  the  precession  rate  4*  is  initially  zero  for  the  case  of  Fig.  11, 
and  4  =  0  <  p.p/2,  the  vehicle  will  precess  in  the  negative  mode  following 
reentry  and  will  remain  in  the  negative  mode  until  the  angle  of  attack 
approaches  trim.  This  is  verified  from  the  computer  result.  Alternatively, 
for  the  case  of  Fig.  12,  4  is  initially  positive  and  greater  than  the  reduced 
frequency  |ip/ 2.  Therefore,  4  continues  to  precess  in  the  positive  mode  and 
remains  positive  as  long  as  the  minimum  value  of  the  oscillation  envelope 
reamins  greater  than  |ip /2 .  The  reduced  frequency  is  also  plotted  in  Fig.  12 
and  remains  less  than  the  minimum  value  of  4  throughout  the  trajectory. 
Consequently,  4  remains  in  the  positive  mode,  as  verified  by  the  computer 
result. 

Figure  13  corresponds  to  a  reentry  in  which  the  vehicle  is  coning 
initially  about  an  axis  inclined  to  the  flight  path.  The  precession  rate  4  with 
respect  to  the  flight  path  is  positive  initially  and  oscillates  with  an  upper  and 
lower  bound.  The  lower  bound  is  plotted  against  the  reduced  roll  rate  in 
Fig.  14  and  is  greater  than  |xp/ 2  initially.  However,  the  roll  rate  increases 
because  of  an  applied  roll  torque,  and  the  reduced  roll  rate  approaches  the 
minimum  value  of  4  after  approximately  8  sec,  which  causes  4  to  become 
unstable.  Simultaneously,  the  lower  bound  of  the  angle-of-attack  oscillation 
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Figure  12.  Motion  History  for  Symmetric  Precession 
About  Flight  Path  at  Reentry 
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Figure  13.  Motion  History  for  Unsyrametric  Precession 
About  Flight  Path  at  Reentry 
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Figure  14.  Precession  Instability 


envelope  approaches  zero  aa  the  upper  bound  of  diverges,  and  4/  reverts  to 
the  more  stable  negative  mode.  The  character  of  the  motion  also  changes 
after  the  instability,  as  shown  in  Fig.  15,  which  is  a  polar  plot  of  the  angle 
of  attack  versus  the  precession  angle  ¥.  Before  the  instability,  the  motion 
is  characterized  by  "flower  petal"  loops  containing  the  origin,  whereas  after 
the  instability,  the  loops  do  not  contain  the  origin  (retrograde  precession). 


Vn.  SUMMARY  AND  CONCLUSIONS 


A  simple  expression  has  been  derived  for  the  angle -of- attack 
convergence  of  a  rolling  reentry  vehicle  subject  to  roll  acceleration  and 
Including  pitch  and  yaw  damping.  The  analytical  result  has  been  shown  to 
compare  favorably  with  computer  solutions  of  the  equations  of  motion.  The 
effect  of  roll  acceleration  on  the  angle-of-attack  convergence  envelope  is 
Identical  to  the  effect  of  pitch  or  yaw  damping  in  the  term  p/p,  which  either 
adds  to  or  subetracts  from  the  damping  coefficient,  depending  on  the  direction 
of  the  roll  acceleration.  The  cumulative  effects  of  damping  and  roll  accelera¬ 
tion  on  the  convergence  envelope  can  be  significant  in  certain  cases.  The 
angle -of -attack  convergence  also  depends  slightly  on  the  mode  of  precession 
motion,  which  depends,  in  turn,  on  the  initial  reentry  conditions. 

It  has  been  shown  that  a  rolling  reentry  vehicle  has  two  modes  of  preces¬ 
sion  motion,  either  of  which  can  persist,  depending  on  the  initial  reentry  con¬ 
ditions  and  on  the  roll  history.  The  positive  mode  io  quasi- stable  and  will 
persist  as  long  as  the  lower  envelope  of  the  precession  rate  oscillation  (40  is 
greater  than  the  reduced  roll  rate  parameter  pp/2.  If  this  parameter 
Increases  because  of  roll  acceleration  until  it  reaches  the  lower  envelope 
of  the  4  oscillation,  a  precession  instability  will  occur,  and  the  precession 
rate  will  change  from  the  positive  mode  to  the  more  stable  negative  mode. 

At  the  point  of  instability,  the  lower  envelope  of  the  angle-of-attack  oscil¬ 
lation  reaches  zero.  Following  the  instability,  the  angle  of  attack  increases 
slightly,  depending  on  the  magnitude  of  the  roll  acceleration  and  on  the  pitch 
a.  '  yaw  damping.  The  character  of  a  polar  plot  of  angle  of  attack  vs  pre¬ 
cession  angle  also  changes  from  "flower  petal"  loops  containing  the  origin, 
before  the  instability,  to  loops  that  do  not  contain  the  origin,  after  the  in¬ 
stability. 
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lit- ASSTRACT 

"^A  simple  expression  is  derived  for  the  influence  of  roll  acceleration  on  the 
angle -of-attack  convergence  of  rolling  reentry  vehicles  with  pitch  or  yaw 
damping.  Also  included  is  an  analysis  of  the  windward -meridian  rotation  rate 
of  the  rolling  vehicle,  which  is  coupled  through  the  gyroscopic  equations  of 
motion  to  the  angle  of  attack  and  roll  rate.  It  is  found  that  two  modes  of  motion 
exist  for  which  the  windward  meridian  is  generally  oscillatory  in  one  mode  and 
rotary  in  the  other  (i.  e. ,  the  vehicle  either  oscillates  about  its  roll  axis  rela¬ 
tive  to  the  wind  or  rotates  continuously  in  one  direction).  The  mode  of  motion 
depends  on  the  initial  reentry  conditions,  and  it  is  shown  that  a  roll-induced 
■  instability  can  occur  whereby  the  motion  changes  from  the  oscillatory  mode  to 
the  more  stable  rotary  mode.  The  analytical  approximations  are  compared 
l  with  computer  solutions  of  the  complete  equations  of  motion,  (y 
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